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Abstract 

The excessive production of gravitinos in the early universe destroys the suc- 
cessful predictions of nucleosynthesis. The thermal generation of gravitinos 
after inflation leads to the bound on the reheating temperature, T^h < 10^ 
GeV. However, it has been recently realized that the non-thermal generation 
of gravitinos in the early universe can be extremely efficient and overcome 
the thermal production by several orders of magnitude, leading to much 
tighter constraints on the reheating temperature. In this paper, we first in- 
vestigate some aspects of the thermal production of gravitinos, taking into 
account that in fact reheating is not instantaneous and inflation is likely 
to be followed by a prolonged stage of coherent oscillations of the infiaton 
field. We then proceed by further investigating the non-thermal generation 
of gravitinos, providing the necessary tools to study this process in a generic 
time-dependent background with any number of superfields. We also present 
the first numerical results regarding the non-thermal generation of gravitinos 
in particular supersymmetric models. 
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1 Introduction 



The overproduction of gravitinos represents a major obstacle in constructing cosmolog- 



ical models based on supergravity |2^. Gravitinos decay very late and - if they are 
copiously produced during the evolution of the early universe - their energetic decay 
products destroy the "^He and D nuclei by photodissociation, thus jeopardizing the suc- 
cessful nucleosynthesis predictions p|, 0- As a consequence, the ratio of the number 
density of gravitinos n^i2 to the entropy density s should be smaller than about 10~^^ 
m for gravitinos with mass of the order of 100 GeV. 

Gravitinos can be produced in the early universe because of thermal scatterings in 
the plasma during the stage of reheating after inflation. To avoid the overproduction 
of gravitinos one has to require that the reheating temperature T^h after inflation is 
not larger than ~ (10^ — 10^) GeV 0. We will come back to this point and present a 
detailed analysis of the thermal generation of gravitinos during reheating. 

However, it has been recently realized that the non-thermal effects occuring right 
after inflation because of the rapid oscillations of the inflaton field(s) provide an extra 
and very efficient source of gravitinos p, p. The helicity ±3/2 part of the gravitino is 
excited only in tiny amounts, as the resulting abundance is always proportional to the 
gravitino mass ^3/2 On the contrary, the helicity ±1/2 part obeys the equation 
of motion of a normal helicity ±1/2 Dirac particle in a background whose frequency is 
a combination of the different mass scales at hand: the rapidly varying superpotential 
mass parameter of the fermionic superpartner of the scalar field whose F-term breaks 
supersymmetry, the Hubble rate and the gravitino mass H. The non-thermal produc- 
tion of helicity ±1/2 gravitinos turns out to be much more efficient than their thermal 
generation during the reheat stage after inflation ^ and it was claimed that the 
ratio for helicity ±1/2 gravitinos in generic supersymmetric models of inflation 

is roughly given by IQ-^Trh/V^/^, where V^/^ ~ 10^^ GeV is the height of the potential 
during inflation. This leads to a very tight upper bound on the reheat temperature, 
Trh ^ 10^(^1/71015 GeV) GeV [|. 

The production of the helicity ±1/2 gravitino has been studied in refs. [^, Q starting 
from the supergravity Lagrangian and in the simplest case in which the energy density 
and the pressure of the universe are dominated by an oscillating scalar fleld $ belonging 
to a single chiral superfleld with minimal kinetic term. An application in the context 
of supersymmetric new inflation models has been recently presented in ref. [^]. 

The equation describing the production of helicity ±1/2 gravitinos in supergravity 
reduces, in the limit in which the amplitude of the oscillating fleld is smaller than the 
Planck scale, to the equation describing the time evolution of the helicity ±1/2 Gold- 
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stino in global supersymmetry. This identification explains why there is no suppression 
by inverse powers of Mp in the final number density of helicity ±1/2 gravitinos and is 
a simple manifestation of the gravitino-Goldstino equivalence theorem: on-shell ampli- 
tudes with external helicity ±1/2 gravitinos are asymptotically equivalent to amplitudes 
with corresponding external Goldstinos for energies much larger than the gravitino mass 
0. This is analogous to the longitudinal W bosons in the standard electroweak model 
behaving as Nambu-Goldstone bosons in the high energy limit. 

This simple observation about the gravitino-Goldstino equivalence becomes crucial 
when the problem of computing the abundance of gravitinos generated by non-thermal 
effects involves more than one chiral superfield. Describing the production of helicity 
±1/2 gravitinos through the equation of motion of the corresponding Goldstino in global 
supersymmetry is expected to provide the correct result in the case in which the scalar 
fields after inflation oscillate with amplitudes and frequencies smaller than the Planckian 
scale. Luckily, this situation is realized in most of the realistic super symmetric models 
of inflation [0. 

The goal of this paper is twofold. In the first part of this work we will still concern 
ourselves with some aspects of the thermal production of gravitinos during the reheating 
stage after inflation. We will perform a detailed analysis of such a process, taking into 
account the fact that reheating is far from being instantaneous. Inflation is followed 
by a prolonged stage of coherent oscillations of the inflaton field. In this regime, the 
inflaton is decaying, but the inflaton energy has not yet been entirely converted into 
radiation. The temperature T rapidly increases to a maximum value and then slowly 
decreases as a~^/^, being a the scale factor of the universe. Only when the decay rate of 
the inflaton becomes of the order of the Hubble rate, the universe enters the radiation- 
dominated phase and one can properly define the reheat temperature T^h- During this 
complicated dynamics, both gravitinos and entropy are continously generated and one 
has to solve a set of Boltzmann equations to compute the final ratio n3/2/s. 

In the second part of this work we will be dealing with the non-thermal production 
of gravitinos during the preheating stage after inflation. 

Our aim is to provide the reader with all the tools necessary to study the helicity 
±1/2 gravitino production in a generic time-dependent background and with a generic 
number of superfields. To achieve this goal, we will derive the master equation of 
motion of the Goldstino in global supersymmetry with a generic number of superfields 
and show that, in the case of one single chiral superfield and amplitudes of the oscillating 
fields smaller than the Planck scale, it exactly reproduces the equation of motion of the 
helicity ±1/2 gravitino found in refs. ^ starting from the supergravity Lagrangian. 
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As a special case, we will concentrate on the case of two cliiral superfields, which is 
particularly relevant when dealing with supersymmetric models of hybrid inflation. We 
will also present the first complete numerical computation of the number density of the 
helicity ±1/2 gravitinos during the stage of preheating for one single chiral superfield. 
This numerical analysis will be performed keeping all the supergravity structure of the 
theory. 

The paper is organized as follows. In sect. 2 we comment about the thermal pro- 
duction of gravitinos. In sect. 3, we show how to derive the equation of motion 
of the Goldstino in a generic time-dependent background, we reproduce the helicity 
±1/2 gravitino equation found in supergravity for one single chiral superfield, we com- 
ment upon the decay rate of the helicity ±1/2 gravitinos and present the numerical 
results regarding the number density of gravitinos in particular supersymmetric models 
containing a single chiral superfield. Finally, in sect. 4 we discuss the non-thermal pro- 
duction of gravitinos for the case of two chiral superfields, which is relevant for realistic 
supersymmetric models of inflation. 

2 Aspects of thermal production of gravitinos dur- 
ing reheating 

At the end of inflation the energy density of the universe is locked up in a combination 
of kinetic energy and potential energy of the inflaton field, with the bulk of the inflaton 
energy density in the zero-momentum mode of the field. Thus, the universe at the end 
of inflation is in a cold, low-entropy state with few degrees of freedom, very much unlike 
the present hot, high-entropy universe. After inflation the frozen inflaton-dominated 
universe must somehow be defrosted and become a high-entropy radiation-dominated 
universe. 

The process by which the inflaton energy density is converted to radiation is known as 
"reheating" ll^]. The reader should rememeber that - even if the process of reheating is 



anticipated by a stage of preheating |T2| - the efficiency of preheating is very sensitive to 



the model and the model parameters. In some models the process is inefficient; in some 
models it is not operative at all. Even if preheating is relatively efficient, it is unlikely 
that it removes all of the energy density of the inflaton fleld. In particular, already 
during the resonant decay of the inflaton fleld, back-reaction processes of rescattering 
|13| always create a sizeable population of inflaton quanta with non-zero momentum 
which do not partecipate in the resonant decay. It is therefore likely that a stage 



during which the inflaton fleld is slowly decaying is necessary to extract the remaining 
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inflaton field energy. This is exactly the stage we are going to analyze in this section. 

The simplest way to envision this process is if the comoving energy density in the 
zero mode of the inflaton (or the soft quanta generated in the process of rescattering 
during preheating) decays into normal particles, which then scatter and thermalize to 
form a thermal background. It is usually assumed that the decay width of this process 
is the same as the decay width of a free inflaton field. 

There are two reasons to suspect that the inflaton decay width might be small. The 
requisite flatness of the inflaton potential suggests a weak coupling of the inflaton field to 
other fields since the potential is renormalized by the inflaton coupling to other fields. 
However, this restriction may be evaded in supersymmetric theories where the non- 
renormalization theorem ensures a cancelation between fields and their superpartners. 
A second and basic reason to suspect weak coupling is that in local supersymmetric 
theories gravitinos are produced during reheating. Unless reheating is delayed, graviti- 
nos will be overproduced, leading to a large undesired entropy production when they 
decay after big-bang nucleosynthesis. 

As we already mentioned, of particular interest is a quantity known as the reheat 
temperature, denoted as Trh- In the oversimplified treatment, the reheat temperature 
is calculated by assuming an instantaneous conversion of the energy density in the 
inflaton field into radiation when the decay width of the inflaton energy, F^, is equal 
to if, the expansion rate of the universe. 

The reheat temperature is calculated quite easily [^. After inflation the inflaton 
field executes coherent oscillations about the minimum of the potential. Averaged 
over several oscillations, the coherent oscillation energy density redshifts as matter: 
Pff) oc a~^, where a is the Robert son- Walker scale factor. If we denote as pi and a/ the 
total inflaton energy density and the scale factor at the initiation of coherent oscillations, 
then the Hubble expansion rate as a function of a is {Mpi = y/SnMp is the Planck mass) 
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Equating H{a) and leads to an expression for a//a. Now if we assume that all 
available coherent energy density is instantaneously converted into radiation at this 
value of a//a, we can deflne the reheat temperature by setting the coherent energy 
density, p^ = pi{ai/a)^, equal to the radiation energy density, pr = {n"^ /?>Q)g^T^^, 
where g^, is the effective number of relativistic degrees of freedom at temperature Tjm- 
The result is 
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T^Mpi =0.2 JT^Mpi. (2) 
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2.1 Thermal production of dangerous relics in the case of in- 
stantaneous reheating 

Under the approximation of instantaneous reheating, the number density of any dan- 
gerous gravitational rehc X is readily solved. The Boltzmann equation reads 

^ + 3Hnx ^ ((Ttotl^^D^ught. (3) 

where atot oc 1/Mp is the total cross section determining the rate of production of the 
gravitational relic and riHght ~ represents the number density of light particles in the 
thermal bath. 

Since thermalization is by hypothesis very fast, the friction term 3Hnx in Eq. (|^) 
can be neglected and using the fact that the universe is radiation-dominated, i.e. H 
t'^ ~ T^/Mp, one finds 

rpi 

The number density at thermalization in units of entropy density reads 

nx , „_2 Trh , . 

— ~ 10 — — . (5) 
s Mp ^ ^ 

As mentioned in the introduction, the slow decay rate of the X-particles is the essential 
source of the cosmological problems because the decay products of the gravitational 
relics will destroy the ^He and D nuclei by photodissociation, and thus successful nucle- 
osynthesis predictions. The most stringent bound comes from the resulting overproduc- 
tion of D + ^He, which would require that the relic abundance is smaller than ~ 10^^^ 
relative to the entropy density at the time of reheating after inflation Q 

^ < 10-12. (6) 

s 

Comparing Eqs. and (^, one may obtain an upper bound on the reheating 
temperature after inflation 

TBj,<(10«-109)GeV. (7) 

If TjiH ~ Mgut, dangerous relics such as gravitinos would be abundant during nucle- 
osynthesis and destroy the good agreement of the theory with observations. However, 
if the reheating temperature satisfies the gravitino bound, it is too low to create super- 



heavy GUT bosons that eventually decay and produce the baryon asymmetry |jT5 . 

In the discussion above, the crucial quantity which determines the abundance of 
dangerous relics after reheating is the reheat temperature Tj^h (or the infiaton decay 



5 



rate through Eq. (^). The reheat temperature is calculated by assuming an instan- 
taneous conversion of the energy density in the inflaton field into radiation when the 
decay width of the infiaton energy is equal to the the expansion rate of the universe. 

However, the reheating process is not instantaneous. Right after infiation the decay 
width of the infiaton is expected to be much smaller than the Hubble rate, ^ H, 
otherwise will violate the gravitino bound. Therefore, the universe undergoes a 
very long period of matter-domination during which the energy density is dominated 
by the oscillations of the infiaton field around the minimum of its potential. These 
oscillations last till the cosmic time becomes of the order of the lifetime of the infiaton 
field. 

In this early-time and prolonged regime of infiaton oscillations, the infiaton is never- 
theless decaying, oc e"'"'**, but the infiaton energy has not yet been entirely converted 
into radiation. The temperature T has the following behaviour. When the infiaton os- 
cillations start and a small portion of the infiaton energy density has been transferred 
to radiation, the temperature rapidly grows to reach a maximum value Tmax and then 
it decreases scaling as a~^/^, which implies that the entropy per comoving volume S is 
created: S oc a^^/^ During this long stage, the universe is not yet radiation- 

dominated. Finally, when t ~ F^^, the infiaton energy density gets converted entirely 
into radiation and the universe enters the radiation-dominated phase. Only at this 
point one can properly define the reheat temperature T^h- Indeed, the reheat tem- 
perature is best regarded as the temperature below which the universe expands as a 
radiation-dominated universe, with the scale factor decreasing as g^^^^T^^, where (7* is 
the number of relativistic degrees of freedom. In this regard it has a limited meaning 
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When studying the production of dangerous relics during reheating, it is necessary 
to take into account the fact that reheating is not instantaneous and that the maximum 
temperature is greater than T^h- This implies that T^h should not be used as the 
maximum temperature obtained in the universe during reheating. The maximum tem- 
perature is, in fact, much larger than Tfj/zQand it is inconsistent to solve the Boltzmann 
equation for the gravitational relics assuming that throughout the period of reheating 
'flight ~ oc and that the reheat temperature Tr/^ is the largest temperature of 
the thermal system after infiation. The goal of the next subsection is to provide a more 
appropriate computation of the number density of dangerous relics generated during 
the process of reheating. For sake of simplicity, we will focus on the gravitino case, but 
our results may be easily extended to other dangerous gravitational relics. 

^As an application of this, particles of mass as large as 2 x 10^ times the reheat temperature Tjih 
may be produced in interesting abundance to serve as dark-matter candidates U%. 
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2.2 A more appropriate approach to thermal generation of 
gravitinos during reheating 

Let us consider a model universe with three components: inflaton field energy, p^, radi- 
ation energy density, pr-, and the number density of the gravitino, We will assume 
that the decay rate of the infiaton field energy density into radiation is T^. We will also 
assume that the light degrees of freedom are in local thermodynamic equilibrium. This 
is by no means guaranteed, but the analysis performed in ref. shows that, even if 



thermalization does not occur, production of gravitinos during reheating is not much 
different . 

With the above assumptions, the Boltzmann equations describing the redshift and 
interchange in the energy density among the different components is 

P0 + 3iJp<^ + r^p^ = 0, 
pR + ^HpR - T^p^ = 0, 

723/2 + SHny2 + (o^tot|?^|) ^3/2 - (%/?) 



0, 



where dot denotes time derivative. Here (dtotl'^l) is the total thermal average of the 
cross section times the M0ller fiux factor giving rise to the gravitino production and 
we have neglected the back-reaction of the gravitino abundance on the radiation energy 
density. The equilibrium energy density for the gravitinos, n^/^, is determined by the 
radiation temperature, T. 

It is useful to introduce two dimensionless constants, a,^ and ax, defined in terms 
of and {cr\v\) as 

= a^M^ {a\v\) = axm'f^ . (9) 

For a reheat temperature much smaller than M^, must be small. From Eq. (^, the 
reheat temperature in terms of and is Trh — cx^J"^ M^Mpi. For = 10^^ 
GeV, must be approximately smaller than 10~^^. 

It is also convenient to work with rescaled quantities that can absorb the effect of 
expansion of the universe. This may be accomplished with the definitions 

$ = p^M^'a"^ ; R = Pru^ ; X = ns/aa^^ . (10) 

It is also convenient to use the scale factor, rather than time, for the independent 
variable, so we define a variable x = aM^. With this choice the system of equations (|^) 
can be written as (prime denotes d/dx) 

$' = -Ci ^ $ 
V^xTR ' 
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The constants Ci, C2, and C3 are given by 
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X^Q is the equihbrium value of X, given in terms of the temperature T as 
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Ml 



(13) 



The temperature depends upon R and g^, the effective number of degrees of freedom in 
the radiation: 

' 30 



T 



1/4 ^1/4 



X 



(14) 



It is straightforward to solve the system of equations in Eq. ([TT|) with initial con- 
ditions at X = xi of R{xi) = X{xi) = and $(x/) = $/. It is convenient to express 
p^{x = xj) in terms of the expansion rate at xj, which leads to 

3 H] 



871 Ml Ml 



X, 



(15) 



Before numerically solving the system of equations, it is useful to consider the early- 
time solution for R. Here, by early time, we mean H ^ T^, i.e., before a significant 
fraction of the comoving coherent energy density is converted to radiation. At early 
times $ $/, and i? ~ X ~ 0, so the equation for R' becomes Thus, 



the early time solution for T is simple to obtain [|17 
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Thus, T has a maximum value of 
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(16) 
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which is obtained at x/xj = (8/3)^/^ 
of TjiH and obtain 



1.48. It is also possible to express 
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Figure 1: The time dependence of the temperature, the radiation energy density, the 
inflaton energy density and the gravitino number density for the case Trh = 10^ GeV 

For an illustration, in the simplest model of chaotic inflation Hj ~ with ~ 
10^3 GeV, which leads to Tmax/Trh ~ 103(200/^,)^/^ for Trr = 10^ GeV. 

For x/xi > 1, in the early-time regime T scales as a~^^^, which implies that entropy is 



created in the early-time regime [1^. So if one is producing gravitinos during reheating 
it is necessary to take into account the fact that the maximum temperature is greater 
than Trh, and that during the early-time evolution, T oc a~^^^. 

The equation of motion of the number density of the gravitino is easily solved nu- 
merically. The results are plotted in Fig. 1. The total cross section for the gravitino 
production is such that ax — 16.6(m3/2/Mp)^ while the total number of relativistic 
degrees of freedom is ^ 230. The inflaton parameters have been chosen to have 
Trh = 10^ GeV, which gives Tmax ^ 10^^ GeV. 

We observe that the quantity (ri3/2/s)/(n3/2/s)finah where s is the entropy density, 
gradually increases with time when is smaller than H, but remains always smaller 
than unity until the inflaton decays at t ~ F^^. This means that most of the gravitinos 
are produced at the last stage of reheating when the inflaton decays and it makes 
sense to talk about Trh- We have also checked that (n3/2/'S)finai approximates well the 
usual estimate one gets neglecting the non-trivial evolution of the temperature of the 
radiation during the period F<^ ^ H. This result can be explained recalling that - during 
the coherent oscillation epoch - the entropy per comoving volume is increasing and the 
abundance of the just-produced gravitinos is continuously diluted by the entropy release. 
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We have also checked that the final number density of gravitinos has a dependence, even 
though weak, on the frequency of the infiaton oscillations M^. This dependence is not 
present in the case of instantaneous reheating, where the number density of gravitinos 
depends only upon the reheating temperature and not on the frequency of the infiaton 
oscillations. 

We conclude that - even though the maximum temperature of 10^^ GeV seems to 
be in contradiction with the usually quoted upper bound of (10^ — 10^°), imposing the 
constraint @ gives the usual upper bound on the reheating temperature Trh ~ 10^ GeV. 
One should keep in mind, however, that the thermal evolution of the universe before the 
epoch t ~ is nonstandard and the physics leading to the bound Tj^j < 10^ GeV is 
much more involved than is usually thought. This observation might be relevant when 
dealing with either a different parameter space for the gravitino, e.g. if the gravitino is 
very light, or with other kinds of dangerous relics. 

3 Non-thermal production of gravitinos and the gravitino- 
Goldstino equivalence 

As shown in refs. p, |^, non-thermal effects occuring right after inflation due to the 
rapid oscillations of the inflaton fleld(s) may lead to copious gravitino production. As 
we noted in the introduction, this occurs because the helicity ±1/2 part of the gravitino 
can be efficiently excited during the evolution of the Universe after inflation. The non- 
thermal generation can be extremely efficient and overcome the thermal production by 
several orders of magnitude, in realistic supersymmetric infiationary models. 

The equation of the helicity ±1/2 gravitino has been found in refs. |^, ^ in the 
case in which the energy density and the pressure of the universe are dominated by 
an oscillating scalar field $ belonging to a single chiral superfield with minimal kinetic 
term. In this section we would like to study the non-thermal production of gravitinos in 
a generic time-dependent gravitational background and for a generic number of chiral 
superfields. 

The equation of the helicity ±1/2 gravitino with a single chiral superfield and min- 
imal kinetic term is identical to the familiar equation for a spin-1/2 fermion in a time- 
varying background with frequency fi, which depends upon all the mass scales appearing 
in the problem, namely the Goldstino mass parameter d%W (where denotes the 

superpotential) , the Hubble rate H and the gravitino mass ^3/2. In the limit in which 
the amplitude of the oscillating scalar field is small, I'I'I ^ Mp, the frequency of the 
oscillations tends to d^yV . The frequency Vt corresponds to the superpotential mass 
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parameter of the Goldstino which is 'eaten' by the gravitino when supersymmetry is 
broken. Therefore, the equation describing the production of hehcity-1/2 gravitinos in 
supergravity reduces, in the hmit of |$| ^ Mp, to the equation describing the time 
evolution of the hehcity-1/2 Goldstino in global supersymmetry and no suppression by 
powers of Mp is present. 

This does not come as a surprise and is in agreement with the gravitino-Goldstino 
equivalence theorem. In spontaneously broken supergravity, the initially massless grav- 



itino acquires a mass through the superhiggs mechanism [[18], [T^, by absorbing the 
Goldstino which disappears from the physical spectrum. Before becoming massive, the 
gravitino, which is a Majorana spin 3/2 particle, posseses only ±3/2 helicity states. 
The Goldstino, a Majorana fermion, provides for its missing (longitudinal) ±1/2 states. 
The equivalence theorem is valid in the limit of large energies compared to ^3/2 where 
the longitudinal component of the gravitino effectively behaves as a spin 1/2 Goldstino 
i- 

Therefore, it appears of advantage to compute the equation of motion of the helic- 
ity ±1/2 gravitino by finding the equation of motion of the corresponding Goldstino 
in global supersymmetry. This procedure is particularly welcome when the problem 
involves more than one chiral superfield and is expected to provide the correct result 
for the number density of helicity ±1/2 gravitinos in the case in which the scalar fields 
after inflation oscillate with amplitudes and frequencies smaller than the Planckian 
scale. This is exactly what is realized in most of the realistic supersymmetric models of 
inflation 



3.1 The equation of motion of the Goldstino in global super- 
symmetry and in a time-dependent background 

Let us now find the equation of motion of the Goldstino when the energy density of the 
background is dominated by a set of scalar fields following the trajectories imposed by 
their equation of motion. We will therefore suppose that the scalar fields are displaced 
from the minima of their potential and are free to oscillate about such minima. This is 
what happens right after inflation and during the preheating stage. 

The identification of the Goldstino requires a generalization of the standard proce- 
dure used in the static case, that is when the scalar fields are sitting at the minima 
of their potential and the cosmological constant vanishes. In the following we will ne- 
glect the expansion of the universe. For the practical purpose of computing the number 
density of helicity ±1/2 gravitinos generated during the preheating stage, this is good 
approximation since the non-thermal production of gravitinos is expected to overcome 
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the thermal generation in those supersymmetric models in which the frequency of the 
oscillations of the scalar fields is much larger than the rate of the expansion of the 
universe and most of the gravitinos are generated within the first few oscillations. Ne- 
glecting the expansion of the universe will also make the identification of the Goldstino 
more transparent. Finally, we will not concern ourselves with a theory charged under 
some gauge group, but suppose that during the evolution of the system some F-term 
is nonvanishing. Our findings can be easily generalized to include the possibility that 
supersymmetry is broken by some (time-dependent) D-term. 
Consider a global supersymmetric theory with Lagrangian 

£ = d^Zid'^z' + '-xa^d^Xi - V{z\ z,) - \ {Wi^XiPLX' + h-c.) , 
V{z\zi)^Wi{W''y. (19) 

Here = H^(z*) is the superpotential, and Xi denote a set the scalar and fermionic 
fields respectively, Zi = {z^Y , Wi = dW/dz"^ and = {1 — 75)/2 is the left-handed 
projection operator. The index i runs from 1 to N, being N the number of multiplets 
and we use the standard convention that the sum is intended when the index i is 
contracted. 

The Lagrangian is invariant under the following set of supersymmetric transforma- 
tions 

5,z' = V2ePLx\ 

5,x' = V2&e, (20) 

where s is the spinor parametrizing the infinitesimal supersymmetric transformation 
and we have defined the matrix 

= [i^z' - {w^y] . (21) 

Given a generic background, supersymmetry is broken if 

{SsX') + 0. (22) 
This happens if the expectation value of the matrix 6* is nonvanishing 

(0*) 7^ 0. (23) 

In particular, for a constant (time-independent) background, one recovers the usual 
condition that supersymmetry is broken if, for some field z*, the F-term is nonvansishing 

(m) 7^ 0. (24) 
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On the other hand, in the case of a time-dependent background, the breakdown of 
supersymmetry comes also from the 



piece in the matrix 6*. This is not surprising since supersymmetry is broken in the 
early universe anytime some form of nonvanishing energy density appears. This is what 
happens during inflation and the subsequent stage of preheating and reheating when 
scalar fields oscillate around the minima of their potential. 

The Goldstone theorem tells us that the Goldstino is easily identified from the 
supersymmetric transformation (pO]) 



n = Qix^■ (26) 

We introduce now the two projection operators 



ete 

el 



fp-Lv. = A -^e- 



where we have defined O^G = ©j^Ofc. The two operators project respectively onto the 
subspace orthogonal to the Goldstone fermion and onto the Goldstone fermion itself. 



Making use of the definition (pif ), we find that, for a background of real fields, 

2 

dt 



^^1 +W 



(28) 



Therefore, the combination G^0 gives the total energy density of the system p which - 
if the expansion of the universe is neglected - remains constant in time. 
The spin 1/2 field Xi can be rewritten as 



9te 

Xt = {P^)vXr (29) 



Xi Xi ~l~ Qj-Q^' 



Notice that the fields xt ^-^e not linearly independent since they satisfy the following 
relation 

QiXi = 0. (30) 

This condition tells us that one of the xt fields may be expressed in terms of the 
remaining (A^ — 1) ones. 
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We now choose the nonvanishing vacuum expectation values of the scalar fields in 
the real direction, Re Zi = (pi/\/2, Im Zi = 0. The equation of motion of the scalar and 
fermionic fields read 



4>i = -W,,Wj, (31) 
h^Xi = kx^ + WijXJ, (32) 

where the dots stand for derivative with respect to time, k = ^ ■ k and we have used 
the plane- wave ansatz Xi ~ e*'^'^ for the space-dependent part. The matrices Bj satisfy 
the following equation 

z/Bi = -W.,jQj. (33) 

Inserting now the decomposition ( |29| ) into Eq. (p^), multiplying by 6j and 0| respec- 
tively and making use of the Eq. (|33D, we get the following equations of motion 

i-f^f] = kg^ri + kx, (34) 
z/^ = -ig^^''r] + k7] + 2elW,,xf, (35) 

where we have defined the following combinations 

G = ^ 

^ ete' 

X = Qlxt (36) 

The matrix Q can be expressed in terms of the energy density p and the pressure p of 
the oscillating scalar fields 

g = — . (37) 

Differentiating Eq. (|3^ ) with respect to time and using Eq. (|35D we find the master 
equation of motion of the Goldstino 

ri + Pr] + i-i^Wi] - 2ik-f^elxi = 0. (38) 

This equation is valid for a generic number of chiral superfields and - because of the 
gravitino-Goldstino equivalence theorem - is expected to provide the necessary tool to 
describe the production of helicity ±1/2 gravitinos during the preheating stage after 
inflation, when the typical energy of the system and field amplitudes (or the frequencies 
of the oscillations of the scalar fields) are large compared to the gravitino mass 171^/2 
and smaller than the Planck scale. We notice the non-trivial result that any time- 
dependent function has disappeared from the k'^t] term; in the ultraviolet regime, Eq. 
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(p^) is solved by plane- waves and particle production shuts off as one would expect from 
general arguments. 

In a static background for which 4>i = 0, we have 6j = 6|, G"'" = 6j = and x = by 
virtue of Eq. (^). The Goldstino equation is solved by plane- waves and - as expected 
- no particle production takes place. 

Let us now consider the special case of one single chiral superfield $ with minimal 
kinetic term. We have = 1, the only physical degree of freedom is the Goldstino and 



Eq. simplifies to 



where 



X^ = X = 0- (39) 



(«7°5o - kg^) V = 0, (40) 



^ = (41) 



ei 

er 

and the matrix 0i is given in Eq. (^Tj) for the case i = 1. Notice that the matrix Q has 
manifestly absolute value equal to unity 

Therefore, it is possible to rewrite g in the following form 

g = e^^^"^, (43) 

By making a field redifinition t] — > exp{i'y^(f)ri, the equation of motion of the Goldstino 
becomes 

i'-f'^rj — kt] — irLcsi] = 0, (44) 

where 

"^'^ = ^ = -8^ 

and we have used the fact that g satisfies the following differential equation 

g 

Eq. ( ^41) is the equation of motion of a spin- 1/2 fermion in a time-dependent background 
given by the oscillating mass rries- We now wish to show that the Eq. ( ^if ) found for 
one single chiral superfield in the limit of global supersymmetry correctly reproduces - 
for amplitudes of the scalar field $ much smaller than the Planck scale - the equation 
of motion of the helicity ±1/2 gravitino found in refs. p, |^ starting from a local 
supersymmetric theory, i.e. supergravity. 



^ = (46) 
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3.2 Non-thermal production of gravitinos in the case of one 
chiral superfield 

Let us first remind the reader of some of the basic results obtained in refs. ^ 
regarding the equation of motion of the hehcity ±1/2 gravitino in the case of one single 
chiral superfield and minimal kinetic term. 

If we start with the supergravity Lagrangian, the single chiral fermion x ^ which 
is the superpartner of the scalar component in the chiral supermultiplet $ - plays the 
role of the Goldstino and can be gauged away to zero, so that no mixing between the 
gravitino ip^ and x is present. Under these circumnstances, the equation of motion of 
the gravitino becomes 

= e^^'P"^,%Vp^„ = 0. (47) 

Here T>p is the covariant derivative and greek letters denote space-time indices. The 
condition V ■ R = gives the following algebraic constraint 

7oV'° = cE7.^', (48) 

i=l 

where the matrix c, in the limit of |$| ^ Mp, reduces to 

c=P±^ = -e. (49) 

p 

Two degrees of freedom may be eliminated using Eq. (^). 

We note that the constraint ( ^8]) may be recovered in the following alternative way. 
The mixing term in the supergravity Lagrangian between the gravitino and the Gold- 
stino is of the form 

xreU,- (50) 

By using the definition ( p6D 1] = Qix the mixing term becomes 

^vQiYQU,. (51) 

Choosing the gauge in which such a term vanishes is equivalent to require that 

Oirel^M = 0. (52) 

This condition gives 'Joip^ = —GJ2^=ilii^\ which coincides with the constraint p8|)p|. 
^The constarint (|5^ ) is easily generalized to the case of many superfields 6i7'^6|'0^i = 0. 
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Because of the antisymmetric properties of the Levi-Civita symbol, the equation 
= does not contain time derivatives and provides another algebraic constraint on 
the gravitino momentum modes. Such a constraint allows to remove two extra degrees 
of freedom and to define two physical Majorana fermion states 1IJ3/2 and which 
may be shown to correspond to the ±3/2 and ±1/2 helicity states respectively, by 
explicitly constructing the helicity projectors in the flat limit The Lagrangian may 
be diagonalized as £ = C3/2 ± ^1/2, where 



3/2 



£1/2 



1p3/2 
^1/2 



5ci 

27°9o ± 2 — 7° - "^3/2a ^3/2, 

h^do + ^— 7° ± ms/2a + kG 
2a 



V'1/2 



(53) 



where a is the scale factor, G = A + i'j^B and A and B are time-dependent functions 



A 



B 



^{^ + ^1/2)' 

A ^ ■ 

- 4— m3/2m3/2 



2m 



3/2 



2— ± 3m3/2^ H — I m 



3/2 



(54) 
(55) 



They may be expressed in terms the pressure and the energy density of the scalar field 
$. Here time is conformal and the line element is ds"^ = a?{T){dT'^ — dx^). 
The diagonal time and space components of the Einstein equation become 



a" 



• 2 

a a 
2— 7 



3M| 



Ml 



dt 



d^ 



dt 



(56) 
(57) 



Using the expression for the gravitino mass m3/2 in terms of the superpotential W 

\wm 



^3/2 = 

we can write the scalar potential V as 



V 



?,\w\- 



m3/2Mp 



am3/2^ 



{5t 



(59) 
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Replacing Eqs. (0) and (|57| ) in Eq. (|59|), one obtains 



m3/, = + - K - 3m3/, + 2- + 6-^3/, (60) 



When this expression for m3/2 is used in Eqs. ( p^ and ([55|), we obtain the ramarkable 
property ||, § 



^2 + 5^ = 1. (61) 



We are now in the position to show the gravitino-Goldstino equivalence explicitly. 
To do so, we neglect the expansion of the universe and the gravitino mass and 
consider the limit |$| ^ Mp. The matrix G has the following limit 0, |^ 

G 1"-^^^ P^^n^ = _gi_ (62) 
P 

The equation of motion of the helicity ±1/2 gravitino therefore reduces to 

(^7°5o - kg^) ^1/2 = 0. (63) 

This equation is exactly reproduced in the global supersymmetric limit by the equation 
of motion of the Goldstino (^). 

One can also use the gravitino-Goldstino equivalence to explain the remarkable prop- 
erty that the matrix G has absolute value equal to unity, by making use of Eq. (|^) . 

3.3 Comments on the gravitino decay 

We bragged about achieving a large number density of helicity ±1/2 gravitinos from 
non-thermal effects and how this phenomenon is strictly related to the fact that Mp 
does not appear in the equation of motions, but then tacitly assumed that the comoving 
number of helicity ±1/2 gravitinos at nucleosynthesis is the same one which may be 
generated during preheating. This issue deserves a closer look because one might think 
that helicity ±1/2 gravitinos promptly decay (or rapidly thermalize), thus not leading 
to a large undesired entropy production when they decay after big-bang nucleosynthesis. 

However, this is not the case; helicity ±1/2 gravitinos do have a decay rate which 
is suppressed by the gravitational coupling Mp^ and is therefore small. This can be 
easily understood in the following way. The helicity ±1/2 components of the gravitino 
field correspond to the Goldstino, which is derivatively coupled to the supercurrent. 
Hence, the total amplitude for the decay rate of the helicity ±1/2 gravitino has to 
be proportional to the mass splitting within the super mult iplets. For example, in the 
present vacuum, where we suppose supersymmetry is broken by some F-term with 
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F ~ m3/2Mp, the coupling between the hehcity ±1/2 gravitino, a fermion / and its 
superpartner / is proportional to {m-— mj) / F. Since {mj—mj) ~ fnl^2^ the coupling 
is suppressed by m3/2/Mp ^. Similarly, the coupling of the helicity ±1/2 gravitino, with 
a gauge boson and a gaugino is proportional to mx/ F ~ Mp^, where mx ~ "^3/2 is 
the gaugino mass. This is the reason why, when dealing with thermal production of 
gravitinos during reheating, the helicity ±1/2 and ±3/2 gravitinos are treated on the 
same ground and have both Mp-suppressed cross sections. 

Right after inflation and during the preheating stage, supersymmetry is badly broken 
by the energy density p stored in the oscillating scalar fields and what measures the 
breaking of supersymmetry is not a simple F-term like in the present vacuum, but the 
parameter V G^B = p^^^. The helicity ±1/2 gravitinos may decay into lighter fermions 
and sfermions through a coupling proportional to Am^/p^/^, where Am^ is the mass- 
splitting in the given light supermultiplet. As supersymmetry breaking is transmitted 
by the gravitational force, at the preheating stage Am? is at most of the order of 
p/M'^ ~ . If the decay is kinematically allowed, the decay rate of the helicity ±1/2 
gravitinos into fermions and sfermions is at most 

p ^ f H \ f VL 



where VL is the (decreasing) time-dependent frequency of the oscillations of the scalar 
fields responsible for the non-thermal production of the helicity ±1/2 gravitinos. Sim- 
ilarly, in the case of decay into gauge bosons plus gauginos, the decay rate is ~ 
{VL/MpfVL < H. 

These estimates are valid as long as the oscillating scalar fields dominate the en- 
ergy density of the universe. The "composition" of the helicity ±1/2 gravitino through 
the Goldstino mixture ( p6D changes with time. During the prolonged stage of coher- 
ent oscillations, the main contribution to the helicity ±1/2 gravitino comes from the 
fermionic superpartners of the coherently oscillating scalar fields and the decay rate ( [6^ ) 
applies. This decay rate is tiny and always smaller than the rate of the expansion of the 
universe; the number density of the helicity ±1/2 gravitinos does not drop during this 
epoch. At later stages, the main contribution to the helicity ±1/2 gravitino is given 
by the fermionic superpartners of the scalar fields whose F-terms break supersymmetry 
in the present vacuum. This means that - when the composition of the helicity ±1/2 
gravitino changes with time - the decay rate will smoothly interpolate between (|6^ ) 



and the more familiar rate F ~ m^i^/Mp. As the decay rate remains smaller than the 
Hubble rate till after the nucleosynthesis epoch, the amount of gravitinos per comoving 
volume generated by non-thermal effects during the preheating stage remains frozen till 
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the age of the universe becomes of the order of Mp/ml^2- At this moment, gravitinos 
decay and their decay products destroy the hght element abundances unless n^/2/s is 
sufficiently small. 



3.4 Numerical results for the case of one chiral superfield 

In this subsection we wish to provide the first complete numerical computation of the 
number density of the helicity ±1/2 gravitinos during the stage of preheating after in- 
flation in the case in which the energy density of the universe is dominated by a single 
oscillating scalar field. It is important to keep in mind that a generic supersymmetric 
infiationary stage dominated by an F-term has the problem that the fiatness of the 
potential is spoiled by supergravity corrections or, in other words, the slow-roll param- 
eter 7] = MpV"/V gets contributions of order unity |Ty]. In simple one chiral field 
models based on superpotentials of the type W = M^$^/2 or = \/A$^/3, super- 
gravity corrections make infiation impossible to start. To construct a model of infiation 
in the context of supergravity, one must either invoke accidental cancellations [BTI], or 



a period of infiation dominated by a D-teim p2[, or some particular properties based 



on string theory pSf. Nevertheless, we are not interested here in the infiationary stage. 



but rather on the subsequent stage of preheating. During this period, it might be that 
the superpotential is well-approximated by a quadratic or cubic expression along the 
oscillating scalar field. 

The equation for the helicity ±1/2 gravitino in the supergravity approach with one 
single chiral superfield has been reduced to a more familiar second-order differental 
equation for a spin-1/2 fermion in a time- varying background in refs. [^, We wish 
to present here a slightly different derivation. Since the matrix G has absolute value 
equal to unity, it is possible to rewrite it in the following form 

G = e'^^"'^, (65) 

where ip is a phase depending upon the conformal time. By making a field redifinition 

'01/2 0'~^^^exp(— i7°(y9)'?/'i/2, the Lagrangian C1/2 simplifies to 

'Ci/2 = i'1/2 - m^sa] i)i/2- (66) 
This is the Lagrangian for a spin-1/2 fermion in a time- varying background with effective 



mass 



if 

^3/2 + [ - 



a 



(67) 
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where (p = —{i/2)^^{G/G) and G = A + i'-f^B is given in eqs. (p^. One can use 
as a guide the recent results obtained in the theory of generation of Dirac fermions 
during and after inflation |2y]. During inflation, since the mass scales present in the 
effective mass rries are approximately constant in time, one does not expect a significant 
production of gravitinos (the number density can be at most 713/2 ~ Hf, where Hj is the 
value of the Hubble rate during inflation). However, in the evolution of the Universe 
subsequent to inflation, a large amount of gravitinos may be produced. During the 
inflaton oscillations, the Fermi distribution function is rapidly saturated up to some 
maximum value of the momentum k, i.e Uk — 1 for k < /^max and it is zero otherwise. 
The resulting number density is therefore Uk ~ ^max- The value of fcmax is expected to 
be roughly of the order of the inverse of the time-scale needed for the change of the 
mass scales of the problem at hand. 

The field ipi/2 can be as usual expanded in terms of Fourier modes of the form 

A/2 = J e-^'-' [uAk,v)ar{k)+Vr{k,r^)bl{-k)] , (68) 

where the summation is over spin and the conditions Vr{k) = Cu^{—k) and ar = 
br are imposed by the fact that the gravitino is a Majorana particle. The canonical 
anticommutation relations imposed upon the creation and annihilation operators may 
be used to normalize the spinors u and v. 

Defining Ur = [u+{'q)il)r{k)^ru-{ri)il)r{k)] and Vr = [rvji.{ri)tpr{k),V-{ri)ipr{k)] , 
where ipr{k) are the two-component eigenvectors of the helicity operators, and using 
a representation where 7° = diag(l,— 1), Eq. (|68|) can be written as two uncoupled 
second-order differential equations for and m_: 

u± + ujI ± i^rricfia)' u± = , (69) 

where, ul = k"^ + m^ga^. In order to calculate the number density, we must first 
diagonalize the Hamiltonian. In the basis of Eq. (|68|) the Hamiltonian is 

H{f]) = J d\Y.{Ek{v) [al{k)ar{k) - br{k)bl{k) 

+ Fk{vM-k)ar{k) + F:{r])al{k)bl{-k)} , (70) 



where the equations of motion can be used to express Ek and Fk in terms of u+ and 

Ek = kRe{u*^u_) + aTTics (l — |w+ 



^Here we choose the momentum k along the third axis and use the representation in which 7'^ 
1 
-1 
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k 

Fk = -{ul - u^_) + amesu+u_ . (71) 

In order to calculate particle production one wants to write the Hamiltonian in terms 
of creation and annihilation operators that are diagonal. To do this one defines a new 
set of creation and annihilation operators, d and 6^, related to the original creation and 
annihilation operators a and through the (time-dependent) Bogolyubov coefficients 
Qffc and Pk, 

d{k) = ak{v)a{k)+f3k{vP{-k) , 

b\k) = -^l[^)a{k) + a*Mb\-k) . (72) 

The Bogolyubov coefficients will be chosen to diagonalize the Hamiltonian. Using the 
fact that the canonical commutation relations imply la^p + 1/3^^ = 1, the choice 
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results in a diagonal Hamiltonian, 

H{'n) = I d\Y.^M [al{k)dr{k) + bl{k)k{k)\ . (74) 

We define the initial vacuum state |0) such that a|0) = b\0) = 0. The initial 
conditions corresponding to the no-particle state are 



u±{0) = \ ; u±{0) = ikuzf{0) =F iamef[U±{0) . (75) 

V CO 

The (quasi) particle number operator J\f = dl{k)dr{k) such that the particle number 
density n is (including the two degrees of freedom from the spin) 

ny2{v) = {0\Af/V\0) = ^^ dkk'lP.f. (76) 

Ti^a-^yrj) Jo 

Let us now consider a quadratic superpotential W = M^^^"^ /2. The supergravity 
potential (|59|) is easily computed for such superpotential. In the limit |$| ^ Mp it 
reduces to = M'^cf? /2, but we have retained its complete supergravity form in the 
numerical analysis. It is useful to write the equation of motions in terms of dimensionless 
variables. We introduce the dimensioneless time f = M^t, as well as the dimensionless 
field X = 

0/00) SO that the scalar field is normalized by the condition Xq = 1. We 
define 0o as the value of the scalar field at the moment when the oscillations start. 

By solving the Einstein Eqs. (0) and (|57D and the equation of motion for the scalar 
field, we have found the time-dependent evolution of nies- It is plotted in Fig. 2 in units 
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Figure 2: The evolution of the parameter mes in units of as a function of time and 
for a quadratic superpotential. The initial condition is (po/Mp = 10~^. 

of and for (po/Mp = 10~^. Notice in particular that at large times, rries tends to 
M^. This is expected since one can verify that, in the limit of |$| <C Mp, rries — —'P/cb 
tends iodlW = M^. 

The result of our numerical integration for the power spectrum of helicity ±1/2 
gravitinos is summarized in Fig. 3 for two different values of initial conditions. Since 
meff changes by an amount ~ M<^ in a time scale ~ M^^, one expects fcmax ~ Af^. 
This expectation is confirmed by our numerical results which indicates a cut-off in the 
spectrum for k ~ M^. 

Finally the ratio of the number density of gravitinos in units of the entropy density 
s is given in Fig. 4 for (po/Mp = 10~^ and in units of the reheat temperature TptH- If 
the mass of the infiaton field is ~ 10^^ GeV as required by the normalization of 
density perturbations, we see that the non-thermal particle production of helicity ±1/2 
gravitinos gives rise to a number density well beneath the bound (^) ^. 

Let us now consider a cubic superpotential W = V^<l>^/3. In the limit |$| <^ Mp the 
potential (|5^ reduces toV = Xcp^/ 4. A special feature of this theory is that the problem 
of gravitino production in an expanding universe can be completely reduced to a similar 
problem in Minkowski space-time by a simple conformal redefinition of the scalar field. 
This explains why the effective mass mefr does not decrease with time, see Fig. 5. 
Furthermore, meg is expected to oscillate with maximum amplitude = v^"\/A0o 

10, P]. This behaviour is well-confirmed by the numerical results given in Fig. 5. 
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Figure 3: The power spectrum of helicity ±1/2 gravitinos for the initial conditions 
0o/-^p = 10~^ and 0o/-^p — 10~^ and for a quadratic superpotential. 
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Figure 4: The ratio rig/s as a function of time for 4>q/Mp = 10 ^ and for a quadratic 
superpotential. 
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Figure 5: The evolution of the parameter rries as a function of time and for a cubic 
superpotential. 
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Figure 7: The ratio Ug/ p^/^ as a function of time for a cubic superpotential. 

The result of our numerical integration for the power spectrum of helicity ±1/2 
gravitinos is summarized in Fig. 6. In this case meff changes by an amount \fX4>o in a 
time scale (-\/A</)o)~^ and one expects k 

max ~ V'A^o- This expectation is again confirmed 
by our numerical results which indicated a cut-off in the spectrum for k ~ a/A0o- 

Finally, the ratio of the number density of gravitinos in units of the entropy density 
p3/4 jg gjyg]^ jj^ pjg^ 7_ Here p indicate the energy density stored in the massless 
oscillating scalar field 0. The result ^3/2 /p^''^ will contradict the bound (|^) by about 
one order of magnitude when the energy density in the scalar field is transferred 
to the energy density of a hot gas of relativistic particles. 

4 Non-thermal production of gravitinos in the case 
of two chiral superfields 

As we already mentioned, constructing a model of inflation in the context of supergravity 
requires some effort. Realistic supersymmetric models of inflation require the mass of 
the inflaton field to be much smaller than the Hubble rate. This is hard to achieve 
in the context of supergravity since supergravity corrections spoil the flatness of the 
inflaton potential [|1^]. However, some exceptions are known and they usually involve 
more than one scalar field. Consider the superpotential 

W = sL^- , (77) 
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where k is a dimensionless coupling of order unity p3, E5|. The canonically-normahzed 



inflaton field is $ = v^lS*!. The superpotential (^) leads to hybrid infiation. Indeed, 
for $ ^ $c = y^sT^/i, = and the potential reduces to V = fi^ plus supergravity 
and logarithmic corrections [^]. If the Kahler potential for the superfield S is minimal, 
the supergravity corections to the mass term of the infiaton field cancel and they do 
not spoil the fiatness of the potential. For $ ^ $c the Universe is trapped in the false 
vacuum and we have slow-roll inflation. The scale /i is fixed to be around 5 x 10^^ GeV 
to reproduce the observed temperature anisotropy. 

When $ = inflation ends because the false vacuum becomes unstable. The fleld 
(j) rapidly oscillates around the minimum of the potential at (0) = 2/i/y^, while the 
fleld $ rapidly oscillates around zero. The time-scale of the oscillations is 0{fi^^). The 
mass scales at the end of inflation change by an amount of order of fi within a time-scale 
~ fi~^. Therefore, one expects /cmax ~ and 113/2 ~ 10~^/i^ 0. After reheating takes 
place, the flnal ratio ^3/2 to the entropy density is 

!^ ^ 10-2^. (78) 

s /X 

This violates the bound in Eq. (|]) by at least four orders of magnitude even if Tr// ~ 10^ 
GeV and imposes a stringent upper bound on the reheating temperature Tptn ^ 10^ GeV 



The estimate ( [78D obtained in ref. [g] was based on the assumption that the results 
on the gravitino production for a single one chiral superfleld model are valid in a theory 
with more than one superfleld. In the following, we wish to show that this assumption is 
justifled. Instead of attacking the problem of the production of helicity ±1/2 gravitinos 
in theories with more than one chiral superfleld from a supergravity point of view, 
we make use of the gravitino-Goldstino equivalence theorem. The identiflcation of the 
helicity ±1/2 gravitino with the Goldstino is well justifled, since the amplitudes of the 
oscillating flelds in the models of supersymmetric hybrid inflation are far below the 
Planck scale. 

We generically denote the two chiral superflelds involved in the generic problem at 
hand by $1 and $2 and by Xi and X2 the corresponding fermionic degrees of freedom. 
The combination x can be expressed, making use of Eq. (|30D , as 



X = elxi + elx2 = -2^7°^X^ (79) 

^2 



where we have deflned A as 



A = WiZ2 - W2Z, = -^7° {&M - &2&\) . (80) 
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Using Eq. (|35|), the master Eq. (|38| ) becomes 

i] + k'^T] + ar] + i-f^k (p^ + a^^"^) r/ = 0, 

where 



A 

Notice that A satisfy the following differential equation 

A (a + a''") 
A " 2 ' 

which is solved by 



a = tl- (e^Gi - elGs) . (82) 



(83) 



A = A(0)e~Jo (84) 

Therefore during the time evolution of the system, A will never vanish. 
Redefining t] exp (— / dt a/2) rj, Eq. (|8lD can be recast in the form 

V + k'v-(^^ + ^y + ^l'k (^te/"^'^*)- = 0. (85) 

Finding an exact solution to Eq. (^), or even studying the problem numerically, goes 
beyond the scope of this paper; we will limit ourselves to outline a standard approxi- 
mation method to estimate the number density of helicity ±1/2 gravitinos. If we define 
Tfj^{t) oc e~*^* to be the solution of the equation at t ^ — oo, i.e. a plane-wave, Eq. ( P5| ) 
can be written as an integral equation 

Vk{t) = rf^{t) + f Vk{t') sin [n{t - t')] Vk{t')dt', (86) 

J — CO 

where 

>4W^(f4)-n»*(c;te/-.)-^. (87) 

Decomposing rik{t) in terms of = [u+(t)ipr{k) , ru_(t)ipr{k)] and vj" = [—ru*_(t)4'r{k), 
u^{t)tpr{k)], where r is the spin index, in the late time region, Eq. ( P^D possesses the 
solution 

ur{t) = aY'u'^it)+Pl^'vlUt), (88) 
where the Bogolyubov coefficient (31/'' is given by 

poo 

= -{i/2n) / e-^'''V;}\t)Vk{t)ur{t) dt, (89) 

J —oo 

and we have let t +oo. If we treat Vk{t) as a perturbative potential, then we can 
solve Eq. ( ^6]) by iteration. To the lowest order in Vk, one has Ur{t) = u™(t) and the 
Bogolyubov coefficient Pl"^' becomes 
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/■oo 

(31' = -(V2fi) / e-2^™t;:?^(t)\4(tX(t) dt. (90) 

J —oo 

The corresponding number of Goldstinos (or, equivalently, helicity ±1/2 gravitinos) in 
a given spin state is therefore 

Nkr = Y.W?- (91) 

r' 

Even though this approximation is expected to offer only part of the information about 
the resonant behaviour of the system, we beheve it provides the right order of magnitude 
for the number density. In typical supersymmetric hybrid models of inflation, like the 



one described by the superpotential (77), the system relaxes to the minimum in a time- 
scale much shorter than the Hubble time ~ Hj^ ^ since the frequency is set by the height 
of the potential V^^^ ^ jj,^ Hi during inflation. The number of particles N^r depends 
upon Vfe(^), the Fourier transform of Vfc(t). Since Vkit) changes by an amount ~ n 
in a timescale ~ Vki^t) rapidly dies out for frequencies VL ^ ^i, Vki^t) oc 1 for 
uj < i^max ~ and zero otherwise. The number of helicity ±1/2 will be ~ ^^^ax ~ l^^i 
confirming the original estimate made in ref. Pj. 
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